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Abstract 

The density of states and the optical conductivity of the doped Hubbard 

model on a Bethe lattice with infinite connectivities have been studied using 

an analytic variant of the Lanczos continued fraction method. The spectral 

weight of the gap states and the position of the chemical potential upon hole 

or electron doping have been studied. We argue that the strong correlation 

effects such as gap states and midinfrared band shown in two dimensions also 

appear in infinite dimensions. 
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The undoped parent compounds of high-T c superconductors are known to have the 
charge-transfer gap between the upper Hubbard band of Cu 3d character and the valence 
band of O 2p character. [|TJ This is mainly due to the strong correlation effects of the 3d 
electrons on the Cu sites. Upon doping, the system undergoes insulator-metal transition 
with creating new states (gap states) around Fermi energy inside the undoped insulator 
gap. H Recent experiments investigating the electronic structure of high-T c superconduc- 
tors to understand the normal state properties have been performed by x-ray photoemission 
spectroscopy (PES) HH] on one hand and x-ray absorption spectroscopy (XAS) f| on the 
other hand. Both clarify the creation of gap states inside the insulator gap upon dop- 
ing. The occurrence of the gap states has also been supported by the mid-infrared (MIR) 
band in the optical conductivity. || There is, however, a controversy in the behavior of 
the chemical potential between the two measurements. [0] PES experiments on hole doped 
La2-xSr x Cu04 L+ $ and electron doped N d2- x Ce x CuO^ || claim that the chemical poten- 
tial remains at nearly the same position, in the gap states even though the doping character 
changes. On the other hand, XAS experiments, which examine unoccupied states, on the 
same materials show two prepeaks which reflect the chemical potential jumps across the 
undoped insulator gap by changing the doping character instead of remaining at the same 
position as shown in PES experiments. This controversy has not been clarified yet. 

In the theoretical side, the simple one-band Hubbard model may be chosen as an effective 
model to describe these features. [|7| However, treating the Hubbard model has proven a stiff 
task because of the difficulty of handling both itinerant and atomic aspects simultaneously 
in a single theoretical scheme. This awkward nature of the Hubbard model makes ordinary 
analytic approach spiritless. For this reason, only numerical approaches have described the 
doped Hubbard model in a relatively successful way. [Q These numerical works performed 
in two dimensions have shown the appearance of the gap states and the shift of the chemical 
potential, which agree with XAS experiments |J . There are other theoretical works showing 
these two phenomena in terms of the different models such as the Emery model and the An- 
derson lattice model. The former has been studied by cluster calculation |§ and the latter 
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by slave boson representation [|K]]. It is believed that the above phenomena probably have 
nothing to do with high-T c superconductivity but with general strong correlation effects. 
Thus we expect that these strong correlation effects could be seen in infinite dimensions, 
too. This leads us to studying the doped Hubbard model in infinite dimensions where an 
analytic approach is possible. The methods developed in infinite dimensions |TT| have ex- 



plained the Mott-Hubbard transition [|T^ -|l4i successfully and even quantitatively described 



the correlated insulator V2O3, recently |[15|jl6| . Large- dimension approaches, however, have 
never been applied extensively to the doped Hubbard model in the ground state showing 
the gap states, the behavior of the chemical potential, and the MIR band in the optical 
conductivity. 

The purpose of this paper is to show that the above mentioned strong correlation effects 
shown in the finite dimensions can also be seen in the infinite dimensions. To do this, we 
obtain the single-particle density of states (DOS) and the optical conductivity for the doped 
Hubbard model on a Bethe lattice with infinite connectivities through a purely analytical 
method. Our analytical results are comparable with the numerical works in two dimensions 
H and XAS experiments [[J. 

The single-particle DOS can be written as W2 



AxM = "4: Jim ReH i: ,(z)| 2= _^ +e , (1) 



where 



E M = —— Ai , (2) 
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where a v = (iLf„, /„)/(/„, /„), A u = (f u , f v -i). The inner product is defined by 

(A, B) = (^ Q \{A, B^}\ty ) where A and B are operators of the Liouville space, B^ is the 
adjoint of B. We choose / = c Jcr and calculate ot u and A u using the recurrence relation 



fv+i — iLf v — a u f u + A u f u _i. \T7l We apply this formalism to the Hubbard model 



H = -t 2 c ^ko + j 53 njanj-tr, (3) 

<jk>cr ja 



where < jk > means nearest neighbor sites and t and U are the hopping integral and the 
on-site Coulomb interaction, respectively. 

Before going to the doped system, we determine the gap width of the undoped insulator 
by getting the single-particle DOS. We suppose the ground state is antiferromagnetic with 
quantum fluctuation 5, i.e., (nj a ) = (1 — 5), (n ka ) = 5, (ni a ) = (1 — 5), and so on. Since a 
substantial amount of quantum fluctuation brings the system in between the Neel state and 
the paramagnetic state, it is hard to treat the intermediate state analytically. We know, 
however, that the antiferromagnetic ordering effect mainly appears in ot u and the band width 
is mainly determined by A„. Thus using the above formulae we get 

a 2 v = -iU(l - 5), a 2 v+i = —iU5 (4) 

for v > and 

rr2r/1 rN 2 . 2U 2 qt 2 5(l - 5) + q 2 t A 

A^U 2 5(l-5) + qt 2 , A 2 = ^_, ) ;^ 

A 2 , + i = U 2 5(l -5)[1 + 0( ^J!— )] ee a, 

A 2u+2 = 2 q t 2 [l + 0( u2 * i _ 6) )} = b, (5) 

for v > 1, where q is the coordination number. To make the kinetic energy finite as q — > oo, 
one usually scales t as -^=, ]ll] and we set = 1 in what follows. The single-particle DOS 



obtained using Eqs. (4) and (5) yields the gap of width 2yU 2 (^ — 5) 2 + (y/a — Vb) 2 which 
corresponds to the charge transfer gap in high-T c cuprates. It is shown in Fig. 1 (a) for 
U = 3 and 5 = 0.2. 

For the doped system, the antiferromagnetic order is quickly broken down upon doping 
and the magnetic phase is not yet clear. As a naive candidate we consider the paramagnetic 
state which is reasonable for the appropriately doped regime. For this state, we may write 
(nia) — ^2 f° r a ^ s ^ e ' w here x is the doping parameter. Then we get 
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<*4„+i = -&1 + (2i/ + l)x + 0(x 3 ) + O(i-)], 
a 4u+2 = - (2i/ + 2)x + 0(x 3 ) + O(^)], 

a 4 ,+3 = + (2z/ + 2)x + 0(x 3 ) + O(^)], (6) 



for v > 0, and 



Al "T (1-x) + 2' A2 -^ (1 _ x2) + 2 ' 

A 3 = T (l - 4x 2 )[l + 0(-)], A 4 = T ±^[l + O(-)], 

A 4 , + i = ^[1 - {(21/ + l)xf + O(x') + 0(±)], 

A 4 , +2 = l + 0(x 2 ) + 0(^), 

A 4 , +3 = ^[1 - {(2u + 2)x} 2 + 0(x 4 ) + 0(^)], 

A 4 , +4 = l + 0(a; 2 ) + 0(^), (7) 

for v > 1. 

The key approximation in getting Eqs. (4)-(7) is the Hartree-Fock approximation which 
corresponds to neglecting the spatial correlations, which is valid in the large-g limit. There- 
fore, in calculating the ground state correlation function, we use a decoupling scheme, e.g., 
(nj a ni a ) = (nj a )(ni a ). The same decoupling scheme has been done on the average of an 
operator containing [if, rij a ] whose average vanishes, since no steady current exists in the 
ground state. 

From Eq. (7), we find a remarkable fact occurring in the doped system. The dimension 
of the Liouville space of the operator Cj a (t) was infinite for the half-filled Hubbard model. 
fT2|j The effective dimensionality for the doped Hubbard model, however, is surprisingly 
reduced. One can see this in Eq. (7). If one neglects higher order terms such as 0(x A ) and 
0(1 /U 2 ) compared with unity, which is valid for small- x and large-?/, the Liouville space 
of Cjcr(t) becomes (2n — l)-dimensional when x = 1/n, where n is an integer except 0. One 
can find vanishing A u at x = 1/n [ITS], then the continued fraction (2) is truncated and 



we get the single-particle DOS as a linear combination of the sharp Lorentzian peaks with 



appropriate weights. The DOS for various hole concentrations for U = 3 are shown in Fig. 
1 (b)-(d). It is clearly seen that the insulator gap is filled with the gap states upon doping 
and the chemical potential is located just below the gap states (low-energy peak) and near 
the top of the valence band of the undoped antiferromagnetic insulator. If we increase hole 
concentration, we see that the low-energy peak increases as much as the reduced spectral 
weight in both lower and upper Hubbard bands. This low-energy peak may be caused by 
the enhanced spin fluctuations at site j due to increased holes around it. This can be seen as 
the first peak of the two prepeaks in XAS experiments. f| For the electron doping, there is 
a reflection symmetry with hole doping because of particle-hole symmetry, so the chemical 
potential lies just above the peak induced by the spin fluctuation and near the bottom of 
the conduction band of the undoped insulator. We show the relation between the position 
of the chemical potential \i and the electronic density (n) for both hole doped, (n) < 1, and 
electron doped, (n) > 1, regimes in Fig. 2. It is clearly shown that the chemical potential 
eventually crosses the edges of the insulator gap when the doping concentration exceeds 
a certain value. Both experimental || and theoretical works PHTTI including present one 
show that the gap states newly appear and the chemical potential shifts across a significant 
portion of the insulator gap when the doping is changed from holes to electrons. 

We now show that the gap states created by doping are responsible for the appearance 
of the MIR band in the optical conductivity. The optical conductivity in infinite dimensions 
is obtained using the following formula [19| . 



<t(u) = <t q J dm' J deA (e)A(e, u')A(e, J + u) f ^ + ^ , (8) 

where A^e) = ^ 2 ~ e2 for the Bethe lattice, /(&/) is the Fermi distribution function and uq = 
nt *2hv ~ wnere a -> N, V are lattice constant, number of lattice sites, volume, respectively, and 
A(e, uj) = — ^ImG(k, uj) = —-Im[u + irj — e — S(c<j)] _1 . Use of the momentum-independence 
of the self-energy in infinite dimensions has been made. We get A(e, uj) using the relations 
between the on-site Green function G(u) = — i5^-(— iu) and the self-energy £(a>) applicable 
to the Bethe lattice, which are £(k>) = to — — f° r paramagnetic state and H a (oS) = 



w _ G-aiyi) _ _ i f Qr antiferromagnetic state. JT3| , P75[] The optical conductivity c(a>) has 
been shown in units of <Jq in Fig. 3. The MIR band is clearly seen and quite similar to the 
experiment || though the doping concentrations are larger than those of experiment. There 



is also a numerical work for the MIR band in two dimensions. 20 



It is quite interesting to compare our analytic results in infinite dimensions with the 
numerical results in two dimensions || and XAS experiments ||. Both theoretical results 
show the following agreements. The doping effect removes some of the spectral weight from 
both the lower and upper Hubbard band and creates new states inside the insulator gap. 
However, the chemical potential lies near the top of the valence and the bottom of the 
conduction band for hole and electron doping, respectively. It agrees with XAS experiment. 
In the optical conductivity, on the other hand, the MIR band is shown in the insulator gap 
upon doping. 

In conclusion, all of the above features obtained analytically in infinite dimensions have 
also shown by the numerical work in two dimensions, which implies that some of the strong 
correlation effects may play a similar role in the physics of both two and infinite dimensions. 
Therefore, we argue that the phenomena treated in this work are the general properties of 
the strongly correlated system described by the doped Hubbard model. 
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Figure Captions 

Fig. 1: The density of states for j- = 3 with hole doping, (a) x — (half-filling); (b) x — y^; 
(c) x = jq] and (d) x — |. The dotted vertical line indicates the position of the chemical 
potential /i. O.lt* has been used as a Lorentzian half- width. 

Fig. 2: Electron density per site, (n) vs /i for ^ = 3. The dashed line is to guide the eyes. 
Fig. 3: The optical conductivities for £ = 3 with different hole dopings. The solid, dash- 
dot, and dashed lines are x — |, |, and with 5 = 0.2, respectively. The arrows indicate 
the MIR bands. O.St* has been used as a Lorentzian half-width. 
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